Bogoliubov speed of sound for a dilute Bose— Einstein condensate in a 3d optical lattice 
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We point out that the velocity of propagation of sound wavepackets in a Bose-Einstein condensate 
filling a three-dimensional cubic optical lattice undergoes a maximum with increasing lattice depth. 
For a realistic choice of parameters, the maximum sound velocity in a lattice condensate can exceed 
the sound velocity in a homogeneous condensate with the same average density by 30%. The 
maximum falls into the superfluid regime, and should be observable under currently achievable 
' laboratory conditions. 
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There have been vigorous activities, and impressive achievements, concerning Bose-Einstein condensates in optical 
lattice potentials recently, both experimentally 0, S H [1 H @ and theoretically |llllllT(ilTllTllTllll Il5|. 
culminating in the observation of the quantum phase transition from a superfluid to a Mott insulator with a condensate 
of 87 Rb atoms 0. 

On the theoretical side, substantial effort has been devoted to calculating the low- lying excitations of a condensate 
in an optical lattice, and thus the speed of sound, within the Bogoliubov theory [lit HH ITsj. In particular, it has 
been predicted that in a one-dimensional optical lattice the velocity of propagation of sound wavepackets decreases 
monotonically with increasing lattice depth [Tgt l2fl |. In this letter, we argue that the situation is different for 
condensates in three-dimensional lattices: In this case the velocity of sound undergoes a pronounced maximum when 
the lattice is made successively deeper, which falls into the superfluid regime, and should be observable under currently 
""O ■ achievable experimental conditions. 

We consider a 3d Bose-Einstein condensate subjected to a simple d-dimensional optical cosine lattice id — 1,3) of 
depth V , 

V(r) = ^-J2cos(2k LXl ), (1) 
J> ■ »=i 

t> : 

where = 2n/X denotes the wavenumber associated with the lattice-generating laser radiation of wavelength A. For 
\Q , d — 1 we study the propagation of sound in the direction of the lattice, while assuming that the condensate remains 
homogeneous in the orthogonal plane. We also assume that the depth Vq of the lattice be sufficiently deep so that 
only the lowest Bloch band needs to be taken into account. 

The velocity of sound propagation in a Bose-Einstein condensate filling the lattice then is given by the standard 
expression (l9ll20| 
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where m* is the effective mass pertaining to the lowest band, and U is the inverse compressibility of the condensate 
O ■ in the lattice. 

^ ' In principle, also the effective mass m* appearing in eq. does depend on the density of the condensate, and 
J> , has to be determined by solving the Gross-Pitaevskii equation for the condensate in the lattice 0, HTl Il3l ITsj . This 
is particularly pertinent if one considers a Id lattice filled with an effectively Id condensate, strongly confined by 
restoring potentials orthogonal to the lattice direction, so that the transversal degrees of freedom are frozen out 
entirely. If then each lattice site carries a comparatively large number of particles, nonlinear effects due to the mean- 
field interaction play a decisive role. In contrast, we focus here on truly 3d condensates, with an average density 
N/V of presently achievable magnitude. Taking the moderate value N/V = 10 13 cm -1 , say, and a lattice constant 
A/2 = 426 nm, one obtains an average occupancy of about 0.8 atoms per unit cell in a 3d lattice. Thus, we define the 
site occupancy N s by writing 

- - (3) 

V (A/2) 3 ' [ ' 

and consider values of N s which are on the order of unity. In this case, the density dependence of m* is negligible, as 
wittnessed by numerical calculations for d = 1 |l9l l20j , so that we obtain quite accurate predictions by invoking the 
effective mass, and the band structure, of the single-particle problem. 
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Moreover, for low site occupancies it is not necessary to solve the Bogoliubov-de Gennes equations in order to 
obtain the elementary excitations, but we may perform the Bogoliubov transformation directly on the basis of the 
single-particle Bloch waves Uk{r). We then find for the inverse compressibility the convenient expression 



U = U U d d r\u (r)\\ (4) 

5 1 I unit coll 



where Q = (X/2) d is the volume of the unit cell, and uo(r) is the Bloch function corresponding to zero quasimomentum, 
normalized such that J d d r |u (r)| 2 = f2, with the integral again extending over one unit cell. The factor 

Anah 2 N 

Uo = TT (5) 

m V 

coincides with the customary interaction energy parameter characterizing a homogeneous Bose-Einstein condensate 
with a density of N particles per volume V, which consists of atoms with (bare) mass m possessing the s-wave 
scattering length a. 

One can intuitively capture the physics expressed by eq. (JIJ by imagining an initially homogeneous 3d condensate 
within which a d-dimensional optical lattice is errected adiabatically; the particles then gradually get concentrated 
around the local minima of the lattice potential. This is described by the fact that the Bloch function uq(t) develops 
a maximum within each well; the integral over the fourth power of this function, taken over one unit cell, increases 
slowly with increasing lattice depth. Obviously, the density reached at the potential minima becomes the higher, the 
larger the lattice dimension d; as a consequence, the compressibility U^ 1 of the gas within the lattice also exhibits a 
significant d-dependence. It is this feature which translates itself into a dependence of the velocity of sound on the 
lattice depth which is markedly d-dependent. 

Before presenting numerical data for the velocity of sound in one- and three-dimensional lattices, it is useful 
to consider the analytically tractable limiting case of a deep optical lattice. Introducing the single-photon recoil 
energy E R = fi, 2 fc£/(2m), and the associated recoil velocity v R = hk^/m, we rewrite eq. © in terms of convenient 
dimensionless ratios, obtaining 



— = I m (6) 
v K V 2E K m* ' 1 ' 

Expanding the lattice potential quadratically around the minima of the cosine wells, the Bloch function Ua(r) can be 
approximated by a superposition of the groundstate wave functions of the corresponding harmonic oscillators, if the 
lattice is sufficiently deep. This harmonic approximation immediately leads to the estimate 

i r / 2 t/ \ ^/^ 

- / d d r K(r)| 4 ~ \j± for V Q /E R » 1 , (7) 

" J unit cell V 4 Eft./ 

where here and in the following the "~"-sign means asymptotic equality. By comparison with exact numerical data, 
we infer that for d = 1 and Vq/E r = 10 this approximation is about 12% too high; the error decreases to 5% when 
V /E R = 25. 

This estimate now allows us to specify the condition for the validity of our approach more precisely: Since the gap 
between the lowest two single-particle Bloch bands amounts to Vq/2, and that gap has to be large in comparison with 
the interaction energy per particle in order to justify the restriction to the lowest band, we require Vq U . Using 
the above approximation in eq. J1J for U, together with the expression for the density, this gives for d = 3 the 
condition 

T T \ 1/4 



»4V2^iV s -. (8) 

Since the ratio a/A of the s-wave scattering length a to the laser wavelength A usually is of the order of 1/100, this 
condition indeed limits the validity of our reasoning to occupancies N s of a few atoms per site, a "site" corresponding 
to the volume (A/2) 3 . 

The theory of the Mathieu equation now allows one to state an approximate expression for the width A of the 
lowest energy band when the cosine lattice is sufficiently deep, namely |2lJ, |23j 



e r ~ yfH \e r ) 



3/4 

exp 



(-2y/V /E R ) . (9) 
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For Vq/E-r = 10, this estimate still is about 18.5% too high; the error becomes less than 5% only if Vq/Er > 90. 
Since a quadratic expansion of the tight-binding cosine energy dispersion relation readily connects the band width 
with the effective mass, 



m* _ 4 1 

we obtain 

-3/4 



(10) 



' Q|) exp(2VWl^) forV /£R»l- (11) 



m 47r 3 / 2 \i?R 

Combining now the deep-lattice estimates and l|ll(l , the velocity of sound © takes the form 



V En \Er 



2 1/4 W^ ^ exp -VW^ (12) 



for d = 1, while for d = 3 we find 



assuming Vq/Eh S> 1 in both cases. The exponential decrease of these velocities with increasing lattice depth is due to 
the increase of the effective mass (|llfl . or, equivalently, the reduction of the tunneling contact between the wells, while 
the different exponents appearing in the prefactors can be traced to the d-dependent inverse compressibilities Q. 

In needs to be kept in mind, however, that when the lattice depth Vq exceeds a critical magnitude, the transition 
from the superfluid to the Mott insulator state occurs Since the Bogoliubov theory treats the interaction 

only approximately, it is incapable of describing large depletions of the condensate, and thus does not incorporate 
the transition 0|. Hence, the Bogoliubov speed of sound (J2J), and the above estimates, are meaningful only in 
the superfluid regime, sufficiently remote from the transition point. Within the tight-binding approximation, and 
assuming unit occupancy N s = 1, in a three-dimensional lattice the transition takes place if the ratio of the on-site 
interaction energy per particle, 

~ = 4rf /■ d 3 rKr) |4 j (w) 



and the hopping matrix element A/4 adopts the critical value 

2X5.8, (15) 




where z = 6 is the number of nearest neighbours of each site. (For high N s , the critical ratio approaches 4:N s z 23].) 
The function w(r) appearing in eq. (|14|) is the Wannier function for the lowest band. Approximating this Wannicr 
function once again by t he g roundstate wave function of the harmonic oscillator corresponding to a quadratically 
approximated cosine well |2J] , we find 



U . (16) 
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Er A \EnJ 

In conjunction with the approximate relation © for the band width, the criterion (|15fl then allows us to estimate the 
critical lattice depth for d = 3 as 

!H k2 K)- 1171 

Actually, this simple estimate appears to be quite reliable: Considering 87 Rb atoms, one has an s wave scattering 
length of a = 6 nm. Choosing A = 852 nm then gives V c /E-r rs 12.3, in perfect agreement with what has been 
observed in the experiment 0]. 

In order to obtain numerical data for the speed of sound which are not restricted to the deep-lattice regime, 
we calculate the exact Bloch waves for the cosine lattice numerically and perform the Bogoliubov transformation, 
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FIG. 1: Velocity of sound for a dilute Bose-Einstein condensate with Uo/En, = 0.036 in a Id optical lattice, as function of 
the lattice depth. The full line is the result of the Bogoliubov theory, evaluated numerically; the dotted line indicates the 
asymptotic prediction 1121 , 

assuming the lowest Bloch state uo(r) to be macroscopically occupied; the speed of sound then is extracted directly 
from the slope of the quasiparticle energies e(k) for low wavenumbers k. Figure^shows the result for d = 1, as function 
of the lattice depth, for Uo/Er = 0.036, as corresponding, according to eq. (JSJ), to the Rubidium data a = 6 nm, 
A = 852 nm, and N s = 1. Comparison with the deep-lattice formula (|12l) shows that this estimate describes the 
numerical data quite well for Vo/Er > 20. Here the speed of sound decreases monotonically with increasing lattice 
depth; the perfect agreement of our numerically calculated curve with the ones obtained before by Kramer et al. |19| , 
and by Menotti et al. poj , clearly underlines the correctness of our reasoning. 

The corresponding data for d = 3 are displayed in fig. El Now the numerical data for comparatively shallow lattices, 
for which the estimate <[T3"| does not apply, show a well-developed maximum at Vq/Er rs 6.1, the speed of sound in 
a lattice of this depth being roughly 30% higher than in a homogeneous condensate with the same average density. 
Since the superfluid/Mott insulator-transition occurs, for the parameters considered, only when Vq/E-h « 12.3, as 
discussed above, this maximum falls well into the superfluid regime, where the simple Bogoliubov approach is still 
reliable. (For d = 2, one finds a rather weak maximum at an even lower Vq/ Er.) 

The appearance of a pronounced maximum of the velocity of sound wavepackcts reflects the competition between 
the slowly decreasing compressibility U^ 1 , and the exponentially increasing effective mass m* , with increasing lattice 
depth. For d = 1, the decrease of compressibility is so weak that the increasing effective mass wins this competition 
right from the outset, resulting in a monotonically decreasing speed of sound. However, for d = 3 the stronger decrease 
of compressibility, resulting from the enhanced concentration of the condensate at the lattice sites, can over-compensate 
the increase of m* at least in shallow lattices, giving rise to a substantial enhancement of the speed of sound before the 
increasing effective mass again diminishes that speed when the lattice is made deeper. An experimental observation 
of this maximum, which should be possible under presently accessible laboratory conditions, would constitute an 
important confirmation of our present understanding of the dynamics of Bose-Einstein condensates in optical lattices. 
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FIG. 2: Velocity of sound for a dilute Bose-Einstein condensate with Uo/Eb, = 0.036 in a 3d optical lattice, as function of 
the lattice depth. The full line is the result of the Bogoliubov theory, evaluated numerically; the dotted line indicates the 
asymptotic prediction 1131 . The Mott insulator state, for which the Bogoliubov theory does not apply, occurs for values of 
Vq/Eb. higher than about 12.3. The maximum of the velocity of sound lies at Vo/-Br ~ 6.1, well in the superfluid regime. 
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